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Abstract 

Based on a recent communication by the present authors the question of 
energy dissipation in magneto hydrodynamical waves in an inflating back- 
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ground in general relativity is examined. It is found that the expanding back- 
ground introduces a sort of dragging force on the propagating wave such that 
unlike the Newtonnian case energy gets dissipated as it progresses. This loss 
in energy having no special relativistic analogue is, however, not mechanical 
in nature as in elastic wave. It is also found that the energy loss is model 
dependent and also depends on the number of dimensions. 
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1. Introduction 

T 1 

Q\ [ Interaction of electromagnetic waves with plasma, especially reflection and attenu- 
ation of electromagnetic waves propagating in plasma layers has attracted a great 
deal of attention in recent years due to their potential applications in reflecting or 
absorbing e.m wave energy, broadcasting signals in microwave frequencies etc pp. 
Moreover, radiation and energy loss in cosmic plasma via Cerenkov radiation [2] 
as also Heavy quark energy loss in a weakly coupled QCD plasma [3] were studied 
in the past. But here we discuss a situation where energy dissipation in the prop- 
^ ■ agation of an electromagnetic wave through a cosmic plasma medium is not due 
to any scattering or friction but due to the expanding background. In two recent 
communications [H [5] we have investigated the propagation of an electromagnetic 
wave through plasma in an expanding background in the framework of higher di- 
mensional spacetime(HD). The present work is the continuation of the earlier two 
and possibly the last one in the series where we specifically study the question of 
energy dissipation of the electromagnetic wave in its interaction with the plasma 
field. While great strides have been made by general relativists to address the issues 
coming out of the recent observations in the field of astrophysics and cosmology 
and despite the fact that more than 90 percent of the cosmic stuff in stellar interior 
and intergalactic spaces is made up of matter in plasma state the much sought-after 
union between the plasma dynamics and general relativity still remains elusive. One 
inhibiting factor against the reunion is possibly the fact that both Einstein's field 
equations and plasma equations are highly nonlinear such that a combination of the 
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two makes exact analytical solutions very difficult to get forcing people the alterna- 
tive route to either numerical analysis or a linearized approximation of the plasma 
equations. However, following the well known (3+1) formulation of general rela- 
tivity by Arnowitt, Deser and Misner(ADM) and its subsequent development and 
applications for a covariant formulation of the equations of magnetohydrodynam- 
ics(MHD)in general relativity by Thorne and Macdonald [6] there has been some 
spurt in activities in general relativistic plasma [TJ, [HI [HI HO] • While classical MHD 
is rather well developed, not a great deal is known about the GRMHD and partly 
because of the usual subtleties in defining physically meaningful frame of reference 
in GR one must be wary of applying the classical results with intense gravitational 
fields. In his electrodynamics of moving bodies Minkowski has given a covariant de- 
composition of the electromagnetic fields, which was later extended by Pham Man 
Quan [11] and most extensively by Lichnerowicz [12] in GRMHD. In a series of arti- 
cles listed earlier some of us studied the interaction of electromagnetic waves with 
plasma, propagation of Alfven wave and also found out the dispersion relations and 
other related properties in an expanding background, for which as a test case we 
have chosen spatially flat Friedmann-Robertson- Walker (FRW) cosmology gener- 
alised to higher dimensions. This is the simplest background, yet it still illustrates 
how the curvature as well as the nonstaticity of spacetime can affect simple MHD 
results. For this simple metric it is possible to split the ordinary electromagnetic 
field tensor, into the ordinary electric and magnetic fields E and B in terms of 
which the field equations are more familiar. For our present work the consequences 
of this simple decomposition is nontrivial. This allows us to make use of the intu- 
ition from the known results of flat MHD to be applied to its general relativistic 
analogue. As the split formalism has been extensively discussed in the literature we 
shall not , for brevity, restate those things here but refer the readers to reference 
[4]. We have attempted the analysis in the framework of HD spacetime which is 
an active area of activity in its attempt to unify gravity with all other forces of 
nature. It also finds increasing applications in brane inspired cosmology as also in 
STM theory [T3]. Moreover, it may also provide an alternative physical explanation 
of the current accelerating era of the universe without bringing in any hypothetical 
quintessencial type of scalar field [H] by hand. For our case it assumes particular 
significance because it is in the realm of early universe that both plasma physics 
and HD are specially relevant. In fact it can be shown that Einstein's equations 
generalized to higher dimensions admit solutions where as the 3D space expands 
with time the extra dimensions shrink in size as to be currently invisible with the 
present day experimental technique. It is further conjectured that some stabilis- 
ing mechanism (quantum gravity may be a possible candidate) stabilises the extra 
space to the planckian scale. So the universe we see today appears manifestly four 
dimensional. 

As mentioned earlier we have in the past rather extensively discussed the propaga- 
tion of electromagnetic waves through different types of perfectly conducting plasma 
medium. For the very simple, conformally metric form chosen for our test case we 
find that most of the wellknown the general relativistic results mimic the analogous 
newtonnian plasma mechanics, except that all the field variables are now non static 
and share the background expansion of the FRW metric. In this work we turn our 
attention to another important aspect of wave propagation i.e., energy dissipation. 



In classical MHD one finds that no power is dissipated during propagation for a per- 
fectly conducting medium. But the situation drastically changes for an expanding 
background where energy dissipation does occur. It is observed that the expanding 
background introduces a sort of dragging force to effect the dissipation, which is not 
to be confused with the usual mechanical drag. In the present work we have calcu- 
lated the energy dissipation using our solutions from the previous works. Although 
we have mainly carried out the exercise in a higher dimensional (HD) background 
we believe that most of our findings, barring some qualitative differences, apply to 
the four dimensions also. 



2. Mathematical formalism 

a. Newtonnian Mechanics 

Before embarking on the question of dissipation of electromagnetic wave through 
a plasma medium in general relativity we try to recapitulate, very briefly, the anal- 
ogous situation in flat space. From any standard textbook on Plasma physics we 
know that in newtonnian mechanics the dielectric constant is hermitian which im- 
plies no damping. This can be easily shown as follows: When propagating through 
a plasma medium the electromagnetic variables, for example, are given by 

E = Re Eoe^-^ (1) 

J = Re J e i(fcia, - Wt) (2) 

Here E and J may depend on lo and k but not on time such that the average 
power dissipation in a cycle is given by 

P = E l J l = l -[ Eo-he 2 ^ + E J* + E*J + E^e' 2 ^} (3) 

where = k.x — ut Again J = crE. Now for homogeneous spacetime the conductiv- 
ity, a behaves like a scalar. But for anisotropic spacetime (say an external magnetic 
field in a particular direction)it is, in general, a tensor of rank two because, as is 
wellknown, the linearized particle velocity 

v = -^(E + -xB ) (4) 

iuum c 

introduces an anisotropic velocity field such that the constitutive relation reduces 
to, Ji = (JijEj and we get for average dissipation over a complete cycle 

<P>=\\ E 0l ^E* 0j + E^jEoj] (5) 

Now, for any arbitrary vectors and matrix A.M.B = B.M T .A where T refers to 
transposition. So the last equation implies that 

<P>= 1 -E* 0l ((a^ + a tJ )E 0j (6) 



We also know from elementary plasma mechanics that the dielectric tensor of the 
plasma medium is related to the conductivity tensor as 

eij = S i:j + (7) 

Now the dielectric tensor for cold plasma (anisotropic in general) is hermitian, as 
can be checked in any standard text book which necessitates that the conductivity 
tensor should be anti hermitian i.e. a a b = So there is no energy dissipation in 

classical plasma dynamics when an electromagnetic wave moves in a plasma medium 
with a static background. 



b. General Relativistic Case 

The situation changes drastically when a similar analysis is carried out in a 
nonflat expanding background of arbitrary dimensions. As mentioned in the intro- 
duction the authors of this report studied [4, 5] the propagation of electromagnetic 
waves in an expanding plasma background in the framework of Einstein's field equa- 
tions both in four and higher dimensions taking spatially flat Friedmann- Robertson 
-Walker metric for simplicity. Using the well known 3+1 decomposition formalism 
of ADM we get the interesting results that the field variables mimic the classical 
special relativistic results except that the sinusoidal vibrations need to be replaced 
by Hankel functions and the field parameters are no longer constants but share the 
background velocity of the embedded metric. Secondly the close resemblance to 
the special relativistic results may be due to the very simple metrics form we have 
considered - the conformally flat FRW line element. With more complicated back- 
ground metric, we suppose, the 3+1 split would yield significantly different results. 
We shall not go into the details of our earlier works here but to make the present 
work more tractable and transparent we need to digress time to time to one of our 
recent papers [4] very briefly and refer to the relevant equations only as and when 
absolutely necessary. 

For our background space we take the (d+1) dimensional generalized FRW space 
time as 

ds 2 = dt 2 - A 2 (dx 2 + dy 2 + dz 2 + dil? n ) (8) 

(n = 5, 6, 7, ... , d ) 

where A = A(t) is the scale function. 

In an earlier work [15] one of us extensively discussed the (d+1) dimensional isotropic 
and homogeneous space time and assuming an equation of state, p = ■yp found the 
scale factor as (p = pressure, p = energy density) 

A ~ (9) 



We then wrote down the Maxwell's equations (for more details see Mcdonald et 



al for 3+1 split) for this metric. Before proceeding further let us ask the perti- 
nent question - why is it not possible to formulate the equations of electrodynamics 
in manifestly covariant form using electric and magnetic four vectors? This is at 
variance with the case of spin in an external electromagnetic field where one can 
define a spin 4-vector s M whose spatial part reduces to s in the proper frame of the 
particle and so subject to the constraint that s^u^ = with the 4-velocity of 
it |16J . A possible answer to this question comes immediately to mind- the electric 
and magnetic fields are not the spatial component of any four vector. It is only a 
very particular combination of their components which form a fully covariant ob- 
ject, the electromagnetic field tensor, F^ v . Only using this tensor the manifestly 
covariant form of the Maxwell's equations can be achieved. However, if one chooses 
a preferred coordinate system it is indeed possible to use the electric and magnetic 
fields as (d+l)-vector (we are here considering a (d+1) dimensional spacetime) and 
then finally write down the Maxwell's equations in a fully covariant form [17]. As 
mentioned earlier the present work investigates plasma physics in curved spacetime. 
To make use of the intuition gained from the conventional plasma dynamics it is 
preferable to split the electromagentic tensor, F^ v into electric and magnetic fields 
E and B in terms of which equations are more familiar. This requires choosing a 
particular set of fiducial with respect to which E and B and other physical quantities 
are measured. In what follows we shall presently see that for our simple background 
metric the electromagnetic field tensor via (3+1) decomposition does decouple as d- 
dim. electric and magnetic field and for the privileged fiducial observers(FIDOs) 
one may write 

= E v u» - E»u v + e^u^Bg (10) 

J» = PeU»+f (11) 



& = (jyu? + A" (12) 

F*^ = B v u» - B»u u + e^E 1 (13) 

where u M is the fiducial (d+l)-velocity and e M „ 7< 5 is the (d+1) dimensional Levi- 
Civita tensor. 

Here the RHS terms are measured by the fiducial observers in the usual manner of 
flat spacetime, and which therefore have the usual physical interpretations and are 
orthogonal to « M whereas the LHS terms are the reconstructed charge- current (n+1)- 
vector, (n+l)-vector potential, $ M etc. 

Moreover the electric current, J M is now the sum of the two terms corresponding to 
the convection current and the conduction current, j M respectively and j^u^ = 
One can invert these relations to get 

p e = -J% (14) 

f = TJ» (15) 
E n = F ^ Uv ( 16 ) 

B» = - l -e^ 5 u v F lS (17) 



= (18) 

A 1 * = Y"u» (19) 

We are now in a position to formulate the general relativistic Maxwell's equations 
for our simple FRW metric to get (see ref.2) 

V.E = Anp e (20) 

V.B = (21) 
BE 

— = KE + cA' l V x B-AuJ (22) 

C/6 

— = KB- cA~ l V x E (23) 
dt 

Qp e 

= Kp e — V.J (charge conservation) (24) 
and finally the particle equation of motion in (d + 1) dimensions as 

DAd ~ lp A^qfE + A-xB) (25) 



Dt * \ c 

or 

^= d -A Kp + q ( E + A l xB ) (26) 
Dt a \ c / 

where 

= + «-V) (27) 
is the convective derivative and the d- momentum 

p = m e Tv (28) 

(m e is the rest mass, T is the boost factor, and v, the d- velocity). 

Here V. and Vx are the ordinary Minkowskian divergence and curl in Cartesian 
co-ordinates. Thus we see that at least for the very simple type of metric chosen the 
electromagnetic field tensor of general relativity is split up and gets decomposed as 
ordinary flat space electric and magnetic field. 

In what follows we shall consider, for simplicity, the small amplitude linear theory 
such that the convective derivative simply reduces to ordinary derivative, 4r. 
We see that the equations(20, 21) have the form familiar from flat- spacetime, 
Lorentz frame electrodynamics. They permit one [18] to charcterise E and B by 
electric and magnetic field lines while the rest have a slightly different form with 
some additional inputs from curved geometry (e.g., A and K terms). 
In this section we investigate the situation where a plasma in thermodynamic equi- 
librium is slightly disturbed through the passage of an electromagnetic wave. We 
assume that an external ambient magnetic field is also present. We, however, as- 
sume the plasma medium to be cold so that the pressure can be neglected when 
considering the particle equation of motion. In stellar systems one often encounters 
situations where relaxation times are much larger than the age of the universe so 
that collisions (hence pressure) may be neglected. The effect of an electric field is not 



generally seriously considered because of the well known Debye shielding effect. The 
general problem of an electromagnetic wave propagating along an arbitrary direction 
with the external magnetic field is given by Appleton and Hartee in the Newtonian 
case when studying the propagation of radio waves in ionosphere. Considering the 
fact that a general solution with arbitrary 9 is very difficult to tackle in an expand- 
ing background with arbitrary number of dimensions we shall restrict ourselves to 
the cases when the electromagnetic wave propagates parallel and perpendicular to 
the magnetic field. However the topic is of great importance in astrophysics and 
space science where electromagnetic wave propagation in magnetized plasma is very 
relevant. 

With the set of equations split to (d + 1) formalism we are now in a position to 
attempt applications in varied plasma phenomena. 

If as usual we set fc^c = u>i (the angular frequency of the wave at some initial time 
t = ti) then we get from the above Maxwell's equations ( see reference 2 for details) 



= ESi i y a + 7) ^ ^^r^r eih . r 

V 7^(1 + 7) -2} 



- w^^ ( -* e ~* &vt ' (29) 

where 

u d = Uit' 1 ^ (30) 

gives a measure of the red shift of the photon due to background expansion. For 
radiation dominated era 7 = 4 , uj<i = u{t d + 1 , so the rate at which the frequency 
decreases is maximum in 4D universe. Moreover damping is greater in radiation era. 

This finding may have nontrivial implications for astrophysics. In an earlier work 
[T9] one of us showed that the process of nucleosynthesis in higher dimensional space 
time is markedly different from that in 4D space time. So like the previous classical 
case the equation (29) may again be written as 

E» = E%(x,t)e i{ - kiX -^ (31) 

with the essential difference that here Eq depends both on space and time and all 
the other physical quantities like k and u share the expansion of the universe and 

Moreover, the exponent (kiX — ut) may be written in a tensorial form as k a x a , where 
the 4- vector k a = (k, —cot). 

Now referring again to our earlier paper we find that for a two component plasma 



2 



v = — E = — E 33 

m e L uji m e L Ud 

The last equation is very similar to the flat space case except that here, uid is not a 
constant but shares the background expansion. With J M = n qv^ we get 

r = J% e i(k - r - ut) (34) 



where 



ji = .^J_M+lr^ (35) 



This relation via equation(24) simplifies to 



ie 2 -2 



J = E (36) 



In a recent communication [4] we have shown that when an electromagnetic wave 
propagates through a plasma medium in a (d+l)dimensional expanding background 
in the presence of an ambient external magnetic field along the z direction the 
dielectric tensor comes out to be 
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e 33 
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Pi (say) 







P2 



(37) 

(38) 
(39) 
(40) 



so that in matrix form 
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(41) 



Here u p is the plasma frequency given by 



<4 = 



b d n e 2 



eB 



-td(l+~/) = 



eB 



(42) 
(43) 



where B, the orthogonal magnitude of the ambient magnetic field is given by B = 

2 

\(B) Z (B) Z \ 1/2 = Bt d ( 1 +^ for our system. Further the constant 



2 d/2 7r d/2 

= (d=i)i! 

2 (d+l)/2 7r (d-l)/2 



(d — even) 
(d — odd) 



(44) 
(45) 



(d-2)!! 

Thus the introduction of the magnetic field generates varied modes transforming the 
dielectric constant scalar e in equation (35) to a second rank tensor e^. Although 



the equations (30) - (40) exactly resemble the analogous expressions in Newtonian 
theory the fact remains that all the frequencies now depend on time rather than 

2(d-l) 

being constant. Further the cyclotron frequency uj c decays as t di - 1 +~''> exactly similar 
to the orthogonal component of the magnetic field. 

A little introspection of the dielectric tensor shows that like the special relativistic 
analogue it is the anisotropy character of the plasma medium which introduces the 
tensorial behaviour of conductivity because if the external magnetic field is switched 
off , P2 vanishes and p\ = p 3 and it becomes a pure scalar. Moreover it is manifestly 
anti hermitian when u is a real. So one would expect that like the previous analysis 
the average energy dissipation will be nil. But it is definitely not the case as the 
following analysis shows. Let us once again recall the energy dissipation expression 

P = E lt J» = \[ E^JS*** + E,^ + E^J* + £ Vo *e" 2i 1 (46) 

where = k.x — ut and E m and Joi are both functions of time due to the expansion 
of the background space time. One may at this stage use the equation (7) to find 
that 

<V = — ~ <W ( 47 ) 

Od 

Following the previous argument and using the usual matrix transformation rules it 
follows that the second and the third terms of the equation (46) give 

\ E* 01 (a* t + a tJ ) = (48) 

This is like the special relativistic case discussed earlier but the similarity just ends 
there. Skipping intermediate steps we shall calculate the average value of the rest two 
terms over a complete cycle in some mathematical details to visualize the differences 
from the previous case. A long but otherwise straight forward calculation now gives 
that 

-2(l+2d) 

< p >= a (d) < t 2Sin2(kx - cut) > (49) 

where aid) = __Es£Ml±2L 
wneie ayu) 7rm e (i+ 7 )<i-2 

The equation(49) decomposes to 

r T -2(l+2d) r T -2(l+2d) 

< p >= a(d) [ Sin 2k t x / * Cos 2utdt - Cos 2^x / t d< - 1+ ^ Sin 2utdt } 

Jo Jo 

(50) 

We have so far analyzed the whole situation in (d+1) dimensional spacetime with 
a general equation of state p = 7p. But, as pointed out earlier, both higher dimen- 
sional cosmology and plasma phenomena are most relevant in the very early phase 
of the universe when the cosmology was in radiation dominated state with 7 = 1/d. 
This input will considerably simplify the already cumberous mathematical expres- 
sions although we believe most of our inferences will be valid in a general equation 
of state also. For the last term in equation(50) we find through equation (24) 



1 _4d_ f -(7 d+1) „ 

— ftd-i I u d-i Sin u du (51) 
Wi J 



where fit^ 1 = u 2 Here f3 = 2m ^ i a d . Also the time period which also is sharing 
the expansion of the cosmos is given by 

T(t) = — (52) 

The equation(51) integrates to 

u d -! / Smu du 7- / W-i Sin — (— (53) 

Jo (a - 1) Jo <*>i VP 

While dealing with the other expression we get almost similar results. The special 
integrals we are dealing here are called Fresnel Integrals which admit the following 
power series expansions that converge for all x 

S(x)= / Sin u 2 du = V (-l) n 7, ^-t tt (54) 

1 Wo ^ ' (4n + 3)(2n+l)! 1 } 

px 00 T 4n+1 

c (X ) = J g«»» a *' = B-i)- (4B+1)(2B) , < 55 > 



The limits of these functions as x goes to 00 are known, being equal to ^f^. Also 
series expansion tells of non zero value of it for any x. Thus we can conclude 
that the average rate of doing work is non zero as an electromagnetic wave passes 
through a plasma medium. This dissipation is unique and has no special relativistic 
analogue coming as it does due to the expansion of the background. As commented 
earlier that with expansion the background the density of the lines of force due to 
the ambient magnetic field gets thinned out which, in turn,results in the apparent 
damping of the magnetic field. It may once again be pointed out that this type 
of damping is not mechanical in nature as one observes in Axionic dissipation or 
that due to friction or viscosity or collision common in natural processes but may 
be termed as expansion inspired damping. 



3. Discussion : 

We have here invoked the wellknown 3+1 split formalism of physics to electrody- 
namics in FRW-like background. Although we have worked out the problem of 
energy dissipation of a plasma wave in HD spacetime we believe that most of the 
findings are valid, at least qualitatively, in the 4D spacetime also. To start with the 
FRW cosmology is chosen for the very simple reason that it is most easy to handle, 
yet it illustrates how the curvature and nonstaticity of the background can change 
the plasma MHD results. In the process we have got the interesting result unlike 
the special relativistic case where no dissipation occurs the curvature and expan- 
sion do introduce a new phenomena of dissipation unique in general relativity only 
without having any classical analogue. This is not caused by the usual mechanical 
friction-type forces as one frequently encounters in classical mechanics. This may 
be interpreted as caused due to the thinning out of the magnetic field lines density 



due to background expansion which results in the attenuation of the magnetic field 
strength and power dissipation. As a future exercise one should consider an inhomo- 
geneous background to check what role inhomogeneity plays in the whole process. 
Generalization to non linear plasma may also introduce interesting physics. 
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